How much ^°Ni can be produced in Core-Collapse Supernovae? : 
Evolution and Explosions of 30 - 100 Mq Stars 



Hideyuki Umeda^ and Ken'ichi Nomoto^'^ 
^Department of Astronomy, School of Science, University of Tokyo, Hongo, Tokyo 113-0033, Japan 
Institute for the Physics and Mathematics of the Universe, University of Tokyo, Kashiwa, Chiba 277-8582, 

Japan 

umedaSastron. s .u-tokyo .ac.jp; nomotoOastron. s .u-tokyo .ac.jp 

accepted for publication in the Astrophysical Journal 



ABSTRACT 

Motivated by the discovery of extremely bright supernovae SNel999as and 2006gy, we have 
investigated how much ^^Ni mass can be synthesized in core-collapse massive supernovae (SNe). 
We calculate the evolution of several very massive stars with initial masses M < lOOM© from 
the main-sequence to the beginning of the Fe-core collapse and simulate their explosions and nu- 
cleosynthesis. In order to avoid complications associated with strong mass-loss, we only consider 
metal-poor stars with initial metallicity Z = Zq/200. However, our results are applicable to 
higher metallicity models with similar C-l-0 core masses. The C-l-0 core mass for the lOOAf© 
model is Mco — 42.6M0 and this is the heaviest model in the literature for which Fe-core col- 
lapse SN is explored. The synthesized ^^Ni mass increases with the increasing explosion energy 
and progenitor mass. For the explosion energy of E51 = i?/10^^erg =30, for example, the '^^Ni 
masses of M(^^Ni) — 2.2, 2.3, 5.0, and 6.6 Mq can be produced for the progenitors with initial 
masses of 30, 50, 80 and 100 Mq (or C-^O core masses Mco= 11-4, 19.3, 34.0 and 42.6 Mq), 
respectively. We find that producing M (^^Ni) ^ 4Mq as seen in SN1999as is possible for Mco > 
34 Mq and E51 > 20. Producing M{^^m) - 13Mo as suggested for SN2006gy requires a too 
large explosion energy for Afco 43M0, but it may be possible with a reasonable explosion 
energy if Mco > 60Mq. 



Subject headings: supernovae: general 
reactions, nucleosynthesis, abundances 

Introduction 



supernovae: individual (SN 1999as, SN 2006gy) — nuclear 



Massive stars with initial masses of M ~ 10 — 
I3OM0 form iron cores at the end of their evo- 
lution, and the collapse of the iron cores trig- 
gers supernova (SN) explosions. During the ex- 
plosion explosive Si-burning produces radioactive 
isotope ^^Ni. The released energy by its decay 
to ^^Fe via ^^Co becomes the dominant energy 
source to power the optical light of most of SNe. 
The ejected mass of the ^""Ni is typically M(^'^Ni) 
=0.07 — 0.15Mq for canonical core-collapse SNe, 
such as SN1987A and SN1993J (e.g., Nomoto et 
al. 2004 for a review). 



The explosion mechanism of core-collapse SNe 
is still uncertain. At present most popular mecha- 
nism is the delayed explosion model. In this model 
the gravitational energy released by the collapse 
is first converted into neutrinos inside a proto- 
neutron star, then the neutrinos heat up the mat- 
ter outside to form an energetic Shockwave. This 
mechanism may work in relatively less massive 
SNe of M < 2OM0 to produce the normal ex- 
plosion energy of ^ 10^^ ergs. 

It has been recognized that there are significant 
numbers of core-collapse SNe, such as SN1998bw, 
2002ap, 2003dh, 20031w, that explode with much 
larger explosion energies than canonical SNe. 
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These energetic SNe are called "hypernovae" , 
whose observational indication is the very broad 
line spectral feature and relatively large ejected 
mass of ^^Ni (e.g., Galama ct al. 1998: Iwamoto 
et al. 1998). Theoretical modeling suggests that 
these SNe are more massive than ~ 20Mq, and the 
compact remnant mass is large enough to exceed 
the maximum neutron-star mass (e.g., Nomoto et 
al. 2004). In modeling hypernovae we have found 
that the synthesized M (^^Ni) increases with ex- 
plosion energies in the core-collapse SNe models 
with a given progenitor mass. However, there 
may be an upper limit for the M('''^Ni), because 
the reasonable range of the explosion energy is 
upper-limited and also larger explosion energy 
yields a larger ^He to ^^Ni ratio in the complete 
Si-burning region (e.g., Nakamura et al. 2001; see 
also Figures 5 & 6 below). 

An interesting question is how much ^^Ni pro- 
duction is theoretically possible for core collapse 
SNe. One motivation to study this question is the 
observability of high-redshift SNe, that are explo- 
sions of massive stars in the early universe. In 
the early universe which has still very low metal- 
licity, it is expected that the initial mass func- 
tion could be top heavy or have a peak around 
100 - 200Mq (e.g., Nakamura & Umemura 
1999; Abel et al. 2002). The stars more mas- 
sive than 14OM0 become the pair-instability SNe 
(PISNe). However, we have shown in our previ- 
ous papers (Umeda & Nomoto 2002, 2003, 2005; 
UN02, UN03, UN05 hereafter) that the nucleosyn- 
thesis patterns of very metal-poor stars do not fit 
to those of PISNe (also, Hcger & Woosley 2002). 
Instead many of the stars correspond to hyper- 
novae or energetic Fe-core collapse SNe. From 
precise comparisons, M (^"^Ni) is estimated to be 
typically 0.1 - 0.4Mq (Tominaga et al. 2007b; 
UN05) and comparable to a typical present day 
hypernovae such as SN1998bw. Such hypernovae 
can be modeled by the explosion of 25 — 5OM0 
stars with explosion energy ~ 10 — 40 x 10^^ ergs. 

In the present days, however, there are very in- 
teresting peculiar SNe, SN1999as (Knop, R.. ct 
al. 1999) and SN2006gy (Foley et al. 2006; Ofek 
et al. 2007; Smith et al. 2007). These SNe are 
extremely bright, and the observational and the- 
oretical preliminary results suggest that M(^^Ni) 
~ 4M0 for SN1999as (Hatano et al. 2001) and 
M(56Ni) > 13Mq for SN2006gy (Smith et al. 



2007; Nomoto et al. 2007). 

Although these estimates are preliminary, these 
M(^^Ni) are much more than those in previously 
known SNe. In Table 1, we summarize some prop- 
erties of the presently known hypernova and their 
candidates; i.e., SN type, estimated progenitor 
mass, C-l-0 core mass and the ejected M(^^Ni) 
obtained by spectral and light curve fitting. 

As suggested by Smith et al. (2007), we have 
tempted to regard those extremely bright SNe as 
PISNe. However, the PISN models are difficult 
to reproduce the broad light curve of SN2006gy 
(Nomoto et al. 2007). In order to explain their 
light curves, much more explosion energy or much 
less ejecta mass is required than the PISNe mod- 
els currently available (UN02; Heger & Woosley 
2002). This provides us a motivation to study 
how much ^^Ni production is possible in the core- 
collapse SNe because hypernova models predict 
much narrower light curve than PISNe models. 

The synthesized ^^Ni mass increases with the 
progenitor mass. In the literature the most 
massive progenitor model available is the M = 
70MQ(MHe = 32Mq,Mco = 28Mq) model in 
Nomoto et al. (1997). Using this progenitor model 
we never obtain M('''^Ni) > 6Mq for the reason- 
able range of the explosion energies. 

Heger et al. (2003) discussed how massive sin- 
gle stars end their life for almost entire mass and 
metallicity ranges. The models relevant to us is 
M < I4OM0 metal poor stars and these are in- 
cluded in their discussion. However, no detailed 
explanation for the progenitor models are given in 
the paper. They implicitly assumed that metal- 
poor M ~ 50 — 30Mq stars, i.e., massive C-l-0 
core, collapse directly to black holes and eject al- 
most no ^^Ni. This may be true for the perfectly 
non-rotaing stars. However, the explosion of ro- 
tating black hole forming stars is uncertain and 
they may even explode as hypernovae with large 
amount of ^^Ni. 

In this paper, we thus calculate several very 
massive metal-poor stellar models up to lOO Mq 
from the main-sequence to the beginning of the 
Fc core collapse. Unlike Heger et al. (2003), we 
explore the possibility that these stars explode en- 
ergetically and using these models we investigate 
the '''*'Ni mass production as a function of the stel- 
lar mass and explosion energy. Such very massive 
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stars experience severe mass-loss when the stel- 
lar metallicities are large. Since the amount of 
mass-loss is very uncertain, we only consider such 
metal-poor models a.s Z = 10^^ = Z0/2OO, for 
which the mass loss rate and thus the uncertainty 
in the stellar mass is relatively small. We note 
that the way of explosion, including the amount 
of the synthesized ^^Ni, is basically independent 
of the existence of He and H-envelopes. Therefore, 
most our results are applicable to higher mctallic- 
ity SN models with the same C-t-0 core progenitor 
masses. 

The explosion mechanism of hypernovae is yet 
unknown. However, the inferred large remnant 
mass suggests that a black hole is formed and 
thus the explosion mechanism is different from 
the delayed neutrino heating model. One pos- 
sible mechanism is the formation of the central 
black hole & accretion-disk system, and the ex- 
traction of the energy from that; this is similar to 
the "micro-quasar" (Paczynski 1998), "collapsar" 
models (MacFadyen & Woosley 1999; MacFadyen, 
Woosley & Heger 2001) and the jet-like explosion 
(Maeda & Nomoto 2003a, b). If this is the case, 
we need to perform 2D or 3D calculations of the 
explosion. For simplicity, however, we perform 
only ID calculations assuming spherical symme- 
try. The spherical explosion model gives the upper 
limits to M(^^Ni) because in the jet-like explosion, 
for example, ^^Ni may be synthesized only along 
the jet directions. 

In section 2, we describe our model and the cal- 
culation method. In section 3 we show the results 
of the presupernova stellar evolutions and the syn- 
thesized Af('''^Ni) during the SN explosions. Sec- 
tion 4 gives conclusions and discussions. 

2. Model and Calculations 

We calculate stellar evolution from the main- 
sequence to the pre- Fe-core-coUapse using the 
same method as adopted in UN02, i.e., a Henyey 
type stellar evolution code which is fully coupled 
to a nuclear reaction network to calculate energy 
generation and nucleosynthesis. Some description 
of the code is also given in Umeda et al. (2000) 
and UN05. 

We adopt the solar chemical composition by 
Anders & Grevesse (1989), the solar He abundance 
Yq = 0.27753, and metallicity Zq = 0.02. We as- 



sume that the initial helium abundance depends 
on metallicity Z as Y{Z) = Yp + {AY/AZ)Z, 
where primordial He abundance Yp = 0.247 and 
AF/AZ=1.5265. In this paper, we only deal with 
Z = 10-4 = .Z0/2OO models. We adopt the em- 
pirical mass loss rate (de Jager et al. 1988) scaled 
with the metallicity as (Z/0.02)°-^ (Kudritzki et 
al. 1989). 

As has been addressed by many authors (e.g., 
Weavcir & Woosley 1993; Nomoto & Hashimoto 
1988; Imbriani et al. 2001) one of the most 
important quantities that affect the pre-SN stel- 
lar structure and nucleosynthesis is the carbon 
abundance after He-burning. This quantity sensi- 
tively depends on the treatment for the core con- 
vection during He burning and the reaction rate 
^^C(a, 7)^^0. Therefore, we take the reaction rate 
as a parameter. According to our previous paper, 
we multiply the ^'^C{a,j)^^0 rate of Caughlan & 
Fowler (1988; CF88 hereafter) by a constant factor 
as shown in Table 2. 

As for the treatment of convection, we assume 
Schwartzshild criterion for convective instability 
and diffusive treatment for the convective mixing. 
The time variation of the abundance Xf of a chem- 
cal element i, is governed by the equation 
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(1) 

where [...]nuc represents the change due to nuclear 
reactions, the second term of the right hand side 
represents the effect of convective mixing. In this 
code we first solve the [...]nuc part simultaneously 
with hydrodynamics, and after convergence cal- 
culate convective mixing by the second term of 
Eq.(l) For the diffusion coefficient D of the con- 
vective mixing we adopt the form derived for semi- 
convection by Spruit (1992): 

D = max(0, V,ad-Vad)fkKt(47rrV)' Q - 3^ /V^, 

(2) 

where (3 is the ratio of gas pressure to total pres- 
sure /3 = Pgas/Pi ^ n is the logarithmic gradient of 
mean molecular weight. 



rfln/i 



" rflnp' 



(3) 
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Kf is the thermal diffusivity defined by 



and fk is the squre root of the ratio of "solute" to 
thermal diffusivity. If molecular diffusivity is used 
to estimate fk, the value becomes very small as 
discussed in Spruit (1992). Because we are not 
certain what are the physical process occurring 
in a star, we treat this quantity as a parameter 
(Saio & Nomoto 1998, unpublished; Iwamoto & 
Saio 1999; Umeda et al. 1999). 

Strictly speaking the above coefficient (Eq.(2)) 
cannot be used for the standard convective region 
where /i gradient is nearly zero. In such a region 
we replace by a small finite constant (say 10~^) 
such that the coeeficicnt D takes nearly the same 
order as the mixing length theory. This formal- 
ism is useful when we study the effects of semi- 
convection and convective over-shootings though 
we do not consider such effects in this work. The 
same coefficient is adopted for our previous works 
since 1999: Iwamoto &: Saio (1999) and their later 
works (e.g., Iwamoto et al. 2004). 

For most of the models in this paper we as- 
sume relatively slow convective mixing {fk = 0.1; 
Umeda et al. 2000). Larger fk leads to stronger 
mixing of nuclear fuel, thus resulting in stronger 
convective shell burning. This effect leads to a 
smaller mass core. However, larger fk leads to 
less carbon fraction after He-burning, which in- 
stead leads to a larger core. Therefore, it is not 
simple to predict the final outcome as a function 
of the parameter fk- 

Once prc-SN progenitor models arc obtained, 
we simulate their SN explosions with a ID piece- 
wise parabolic method (PPM) code by putting ex- 
plosion energy into some mesh points just below 
the mass cut. Here we do not specify the explo- 
sion mechanism and assume spherical symmetry 
for simplicity because the hypernova mechanism 
is highly uncertain. The detailed nucleosynthesis 
during the explosion is calculated by postprocess- 
ing as described in UN02. 

3. Results 

3.1. Pre-Supernova Evolution 

In Table 2, we summarize the parameters and 
some properties of our pre-SN progenitor models. 



The first three columns arc the initial mass, con- 
vection parameter fk and the factor multiplied to 
the CF88 ^^0(0;, 'y^^O rate. The next six columns 
are the pre-SN progenitor final-mass with mass- 
loss (M/), He core mass, C-l-0 core mass, C/0 
mass fraction after He burning, Fe core mass and 
the binding energy -Bbin of the layer above the 
Fe-core, respectively. Here the He, C-t-0 and Fe 
core masses are defined by the mass coordinate 
where the fractions of hydrogen and helium are 
X(H)< 10-3, X(He)< 10-3 and electron mole 
fraction Yg < 0.49 in the progenitor models, re- 
spectively. 

The total amount of mass-loss is larger for more 
massive stars. However, with this low metallic- 
ity the total lost mass is not so largo even for 
M = IOOM0 and a large fraction of the hydro- 
gen envelope still remains. Therefore, such a star 
becomes SN II unless it is in a binary system with 
its envelope being stripped off by binary interac- 
tion. 

In the models shown in this paper, the Fe core 
mass is larger for more massive stars. This is 
usually the case but not always, because the Fe 
core mass depends on the C/0 ratio after He- 
burning. For a larger C/0 ratio, convective C- 
burning is stronger, thus forming a smaller Fe- 
core (e.g., Nomoto & Hashimoto 1988). The C/0 
ratio is generally smaller for more massive stars, 
because the ^'^C{a,j)^^0 reaction is more effi- 
cient than the 3a reactions at higher temperatures 
and lower densities. However, the C/0 ratio also 
strongly depends on the convective mixing at the 
end of He-burning. 

The evolutions of M < IOOM0 stars are all 
similar in the sense that all these stars form Fe- 
cores to initiate Fe-core collapse. Nevertheless, we 
find a difference in the evolution of the stars more 
massive than 80M©. These stars experience "pul- 
sations" at the end of Si-burning (see also Heger 
et al. 2003). 

In Figure 1 we show the evolution of the cen- 
tral density and temperature for the 30 and 9OM0 
models. It is well known that a more massive star 
has a larger specific entropy at the center. Thus 
the temperature of the more massive star is higher 
for the same density. Along with this fact we find 
that the central temperature and density evolution 
is very different between the 30 and 9OM0 models 
during the "Si" -burning stage, which takes place 
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around the temperature Tg — 2.5 — 4. Here "Si" 
includes not only ^®Si but also all the elements 
heavier than Si and lighter than Fe. The cen- 
tral temperature and density of the 90Mq model 
oscillate several times. This is because in such 
a massive star radiation pressure is so dominant 
that the adiabatic index of the equation of state 
is close to r ~ 4/3, and thus the inner core of 
the stars easily expands with the input of the nu- 
clear energy released by "Si" -burning. Once it 
expands, the temperature drops suddenly, central 
"Si" -burning stops and the stellar core turns into 
contraction. Since only a small amount of "Si" is 
burnt in each cycle, these pulsations occur many 
times. Note that as shown in Figure 1 the density 
- temperature track is very close to (but outside 
of) the electron-positron pair-instability region" 
where T < 4/3 (e.g., Barkat et al. 1967). 

A similar type of oscillations has been reported 
for metal-poor stars of IOOMq >M> UOMq in 
Heger et al. (2003) as "Pulsational PISNe". They 
discussed that these stars encounter pair instabil- 
ity, leading to violent pulsations, but not complete 
disruptions. Each pulse can have up to several 
2^q5i ergs, and only outer layers of the star are ex- 
pelled. In our model for the 80 — IQQMq stars we 
have not found such mass-loss probably because 
the mass-range is lower and also we have attached 
H-envelope that absorbs and relieves the pulsa- 
tions. 

We have found from the study of 80— lOOM© 
stars that the number of the oscillations depends 
on the convectivc parameter fk'. larger fk in- 
creases the number of the oscillation. This is be- 
cause for a larger fk fresh Si is moiQ efficiently 
mixed into the central region, which increase the 
lifetime of this stage. The amplitude of the tem- 
perature and density variation is larger for more 
massive stars, which suggests that more and more 
drastic oscillations could occur for stars of M > 
1OOM0. During this stage, the specific entropy 
at the center is reduced (or density is increased) 
and the Fe-core, which is defined as a region with 
Ye < 0.49. grows. However, the pulsations does 
not much affect on M{^^m). 

3.2. Pre-supernova Density Structure 

As described in the next subsection, the amount 
of produced ^^Ni is closely related to the den- 
sity structure, more specifically the enclosed mass 



(Mr) - radius (r) relation, of the pre-SN progen- 
itors, where Mr is the mass contained within a 
sphere of radius r. Figure 2 shows this Mr - r re- 
lation for our progenitor models. For reference, in 
Figure 3 and 4 we show distributions of density p 
and electron mole fraction Y^, respectively. We see 
that the Mr - r curve is steeper for more massive 
stars, which means more massive stars are more 
centrally concentrated. In some cases, however, 
less massive stars have a steeper Mr - r relation 
if the C/0 ratio after the He burning is smaller, 
which leads to a larger core mass. 

3.3. Explosive Nucleosynthesis of ''^'Ni 

The products of explosive burning are largely 
determined by the maximum temperature behind 
the shock, Tg. Material which experiences T^ g = 
Tg/lO^ K > 5 undergoes complete Si-burning, 
forming predominantly ^^Ni. The region after the 
shock passage is radiation dominant, so that we 
can estimate the radius of the sphere in which ^^Ni 
is dominantly produced as 

Rm 3700£;*5(^ km, (5) 

which is derived from E* = 4nRf^^aT^/3 with 
Ts^g = 5 (e.g., Thielemann et al. 1996). Note that 
the deposited energy E* is approximately equal 
to the sum of the explosion energy E'exp and the 
binding energy {Eun > 0) of the progenitor, i.e., 

E* ~ i^exp + Ebin- 

We show in Figure 2 the radii which correspond 

to E^i = 1, 10, 50 and 100, where the solid circles 
on the Mr - r curves indicates the outer edge of 
the core, i.e., the masses of Fe-cores. This figure 
suggests that the ^^Ni synthesis is larger for larger 
E* and for the models with steeper Mr - r curves, 
i.e., more massive stars. The mass between the 
outer edge of the Fe-core and r = i?Ni gives the 
crude upper-limit to M(^^Ni) in the ejecta. This is 
an upper-limit because the mass-cut, that divides 
the compact remnant and ejecta, can be far above 
the Fe-core, and also the mass fraction of ^^Ni in 
the complete Si-burning region can be significantly 
smaller than 1.0 especially for the very energetic 
explosion (see Figure 5 & 6 and discussion below). 

In order to obtain the synthesized M(^^Ni) 
more precisely, we perform explosion simulation 
using the ID PPM code and detailed postprocess 
nucleosynthesis calculations (see e.g., UN05). Fig- 
ures 5 & 6 show examples of the results obtained in 
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this way. These are the post-explosion abundance 
distributions in the SOM© and 100 M© models 
for the explosion energy of E51 = 30 (left) and 
E51 = 100 (right), respectively. As shown in these 
figure, for a more energetic explosion, a— rich free- 
zout is stronger so that the mass fraction ratio of 
*He/^^Ni increases in the complete Si-burning re- 
gion. 

The results are summarized in Table 3. Here, 
the models A-F are the same as in Table 2. For 
each model, the results with various explosion en- 
ergies are shown. The upper limit of the ^®Ni mass 
is shown as ^^Niup in the last column. This is the 
upper limit, because the mass-cut is assumed to 
be equal to the Fe-core mass given in Table 2. In 
reality, though it depends on the model assump- 
tionfl^ the mass-cut is likely larger than the Fe-core 
mass, then the ejected ^^Ni mass could be smaller. 
^^Niup is also summarized in Figure 7. 

In Table 3 we show the compact remnant mass 
Mrem for reference. Mi-cm corresponds to the mass- 
cut when 0.07 Mq of '^^Ni is ejected. If M{^'^m) 
is larger, therefore, the remnant mass would be 
smaller than M^cm shown in Table 3. 

In the 7 and 8th columns in Table 3, we also 
show M(^^Ni) (almost exactly the same as the 
ejected Fe mass for these low metalicity models) 
when the mass ratios of Mg/Fe and 0/Fe are the 
solar ratios, i.e., when [Mg/Fe] = and [0/Fe] 
= 0, respectively. We show these quantities, be- 
cause in the observations of metal-poor halo stars, 
most stars have the abundance of [Mg/Fe] > 
and [0/Fe] > 0. This means that if SNe yield 
[Mg/Fe] < or [0/Fe] < 0, or producing larger 
M (^^Ni) than in the 7 and 8th columns, such SNe 
should not be dominant in the early Galaxy. We 
note that if we constrain the ejecta abundances 
as [0/Fe] > 0, ^^Ni ejection of more than 3M© is 
not possible for all the models considered in this 
paper. 

We, however, should make a caution that the O 
and Mg yields from massive stars depend on the 
treatment of convection (e.g., fk) and the adopted 



^If the explosion is strictly spherical and the location of the 
energy deposition is uniquely determined, we can uniquely 
determine the mass-cut as a function of -Boxp- However, in 
an aspherical model or in a jet-like explosion with energy 
injection by continuous jets (e.g., Tominaga et al. 2007a), 
the "effective mass-cut" may not simply the function of 



"'^^C(q!, 7)"'^^0 reaction rate. The model depen- 
dency of O yield would be small, at most a factor 
of 2, but that of Mg yield would be much larger 
and the yield can be a factor of 5 different at most. 
Generally the Mg yield increases when the C/0 
ratio after the He-burning is larger because Mg is 
the C-burning product. 

We find from Table 3 that if relatively less mas- 
sive SNe (M < 30Mo) produce ^^Ni as much as 
2 Mq, the ejecta is too Fe-rich to be compatible 
with metal-poor star abundance. On the other 
hand, the ejection of Af(^®Ni) ~ 2Mq from such 
massive stars as M > 8OM0 does not contradict 
with the abundance of metal-poor stars. 

We note that large M(^'^Ni) (> 4Mq) mod- 
els are quite energetic, so the Zn/Fe ratios in the 
ejecta are relatively large, exceeding the solar ra- 
tio. This is a distinctive feature from any PISNe 
models which do not produce large Zn/Fe (UN02). 
In the energetic explosions Co/Fe ratios are also 
relatively high. However, as shown in UN05, the 
Co/Fe ratio is not much enhanced for Yg < 0.5 
even for the energetic explosions. In the present 
models, indeed, the Co/Fe ratios are less than the 
solar value. More detailed study about nucleosyn- 
thesis and its parameter dependences are given 
elsewhere. 

4. Conclusions and Discussions 

Motivated by the discovery of unusually bright 
SNe 2006gy and 1999as, we have investigated how 
much Af(^^Ni) can be synthesized by core-collapse 
massive SNe with the initial mass of M < lOOM©. 
Observed properties of SN1999as and 2006gy sug- 
gest that the ejected ^^Ni masses are A/(^^Ni) 
~ 4Mq (Hatano et al. 2001) and 13 Mq (Nomoto 
et al. 2007). Previously the only known SNe 
model that can produce such a large M(^^Ni) is 
the PISNe which are the explosion of 140-300Mq 
stars. However, the light curves of PISNe may be 
too broad to reproduce the observed light curves 
of SNel999as and 2006gy, suggesting that a less 
massive star explosion (Nomoto et al. 2007). 

We, therefore, have calculated the evolution of 
several very massive stars with the initial masses 
of M < lOOAf0 from the main-sequence to just 
before the Fe-core collapse. In order to minimize 
the uncertainty in the mass-loss, we have adopted 
low-metal {Z = Z0/2OO) models throughout this 
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paper. However, most results, especially the esti- 
mate of the ejected masses of ^^Ni and ^^O, are ap- 
plicable to higher metallicity models with similar 
C+O core masses. Using these newly calculated 
progenitor models, we simulated the SN explosion 
of such massive stars and calculated the mass of 
the synthesized '"'^Ni, and other elements. 

We find that the M > 80Mq stars "pulsate" 
during the central Si-burning stages. The exis- 
tence of this stage makes the lifetime of the stars 
longer and makes the Fe-core grow. However, the 
Fe-core growth during this stage is not the main 
reason why these massive stars can produce such 
large M(^^Ni). In this sense the existence of this 
stage is not critical for the production of the large 

We have shown that the synthesized M(^^Ni) 
mass increases with the increasing explosion en- 
ergy, Eexp, and the progenitor mass M. Among 
them the effect of M is more important than the 
energy because more massive stars tend to have 
steeper Mr - r curves in the density structure of 
the progenitors, so that M(^^Ni) increases steeply 
with M for the same £^exp- Larger i?exp leads to 
a larger ^^Ni producing region according to the 
Equation (1); however, this also leads to a smaller 
Ni/He ratio in that region. As a result, a large 
i^exp does not result in so large M(^^Ni). For E^i 
= 30, M(56Ni) = 2.2, 2.3, 5.0, and 6.6 Mq can 
be produced for M = 30, 50, 80 and 100 Mq, re- 
spectively, or for C-l-0 core masses of Mco= 11-4, 
19.3, 34.0 and 42.6 Mq, respectively. 

It is possible to produce M(^'^Ni) ~ 4M0 as 
seen in SN1999as if Mco ~ MMq and £^51 > 20. 
On the other hand, producing M{^^Ni) > ISM© 
as seen in SN2006gy is challenging for core collapse 
SN models if Mco < 43M0. In the present study, 
the most massive model has Mqo^ 42.6 Mq. For 
this model, the explosion energy of E51 ~ 200 is 
required to produce M(^^Ni) ~ ISM© but such 
i?5i may be unrealistically large, although there is 
no observational constraints. 

Nevertheless, we should note that if we simply 
extrapolate our results to more massive stars up 
to M ^ ISOMq, Mco is as large as 60Mq. Such 
a model would undergo core-collapse and produce 
M(^^Ni) ~ I3M0 without too large explosion en- 
ergy. Currently we are constructing such an ex- 
treme model, and will present the results in a forth 



coming paper. We are also studying very massive 
stars above 3QQMq for which much study has not 
been done yet. Those stars may be too heavy to 
explode as PISNe and form blackholcs (e.g., Fryer 
et al. 2001). The collapse of these stars may 
or may not accompany aspherical mass ejection 
and potentially eject - IOMq '"^^Ni (Ohkubo et al. 
2006), though the stellar metallicities must be very 
low in order to leave sufficiently massive cores. 
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tific Research (18104003, 18540231) and the 21st 
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Table 1 

Properties of hypernovae and their candidates. 



Name 


Type 


^51 


Mi„i/M0 


Mc+o/Mq 


M(^''Ni)/M0 


Reference 




SN1997of 


Ic 


20 


~ 35 


10.0 


0.15 


Iwamoto ct al. 


2000 


SN1998bw 


Ic 


r-^ 40 


~ 30 


13.8 


0.5 


Iwamoto ct al. 


1998 


SN2002ap 


Ic 


^ 7 


~ 21 


3.3 


0.10 


Mazzali ct al. 


2002 


SN2003dh 


Ic 


~ 30 - 50 


~ 35 - 40 


4.5 


0.35 


Mazzali ct al. 


2003 


SN1999as 


Ic 


~ 20 - 50 


~ 60 - 80 


~ 20 - 30 


4 


Hatano ct al. 


2001 


SN2006gy 


Ib/c ? 


? 


7 


7 


> 13 


Nomoto ct al. 


2007 


SN1987A 


II 


1 


--18-20 


3 


0.07 


Nomoto et al. 


2004 



Note. — Supernova name, type, expected explosion energies in 10^^ erg, initial progenitor masses, C+O 
core masses, ejected ^®Ni masses, and the references are shown. One canonical SN (1987A) is also shown 
for comparison. 
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Table 2 

Parameters and properties of the pre-SN models. 



Model 


M 


fk 


CF88X 


Mf 


A^(Ho) 


M(CO) 


C/O 


M(Fe) 






20 


0.05 


1.4 


19.94 


6.76 


4.39 


0.19/0.80 


1.51 


1.63 




20 


0.10 


1.7 


19.93 


6.45 


5.29 


0.16/0.83 


1.55 


1.38 


A 


25 


0.10 


1.0 


24.87 


8.35 


6.33 


0.29/0.70 


1.52 


1.83 


B 


30 


0.10 


1.0 


29.61 


13.1 


11.4 


0.30/0.68 


1.86 


1.87 


C 


50 


0.50 


1.0 


48.55 


21.8 


19.3 


0.16/0.79 


2.21 


3.67 


D 


80 


0.10 


1.0 


75.81 


36.4 


34.0 


0.20/0.74 


2.75 


5.46 


E 


90 


0.50 


1.3 


83.53 


37.4 


36.5 


0.087/0.82 


3.15 


4.73 


F 


100 


0.10 


1.0 


83.14 


46.5 


42.6 


0.23/0.72 


3.22 


7.19 



Note. — Each column is the initial mass (M), convection parameter (/fc), the factor multiplied to 
the CF88 ^^C(a,7)^^0 rate, pre-SN progenitor mass with mass-loss (M/), He core mass, C+O core 
mass, C/O mass fraction after He burning, Fe-core mass and the binding energy Ebin- Here, all 
masses are shown in the units of Mq and the energy is in 10^^ erg. These models are all for the 
metallicity Z = Zo/200. 
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Table 3 

The masses of Mg, O and upper mass limits of the ^^Ni in the ejecta of core-collapse SNe 

AS A function of PROGENITOR MASS AND EXPLOSION ENERGY. 



Model 


M 


^51 




Mg 


O 


^*Ni[Mg/Fc] = 


''''Ni[o/Fo]=o 


'"'Ni„p 


A 


25 


1 


2.11 


0.14 


3.50 


0.34 


0.46 


0.56 


B 


30 


1 


2.60 


0.41 


5.45 


0.98 


0.72 


0.69 






10 


3.57 


0.41 


4.68 


0.98 


0.62 


1.43 






20 


4.25 


0.36 


4.27 


0.86 


0.56 


1.83 






30 


4.77 


0.31 


3.86 


0.75 


0.51 


2.19 






50 


5.60 


0.26 


3.36 


0.62 


0.44 


2.52 


C 


50 


1 


2.46 


0.79 


12.52 


1.90 


1.65 


0.83 






10 


4.22 


0.73 


11.67 


1.77 


1.54 


1.48 






30 


5.44 


0.79 


10.73 


1.91 


1.42 


2.28 






50 


6.26 


0.80 


9.96 


1.94 


1.32 


2.74 






70 


6.90 


0.80 


9.49 


1.92 


1.25 


2.97 






100 


7.62 


0.78 


8.91 


1.87 


1.18 


3.19 






150 


8.62 


0.76 


8.16 


1.82 


1.08 


3.60 


D 


80 


1 


6.08 


1.41 


18.71 


3.40 


2.47 


2.14 






10 


8.00 


1.35 


17.45 


3.25 


2.30 


3.36 






30 


10.61 


1.28 


15.31 


3.09 


2.02 


4.99 






50 


11.89 


1.20 


14.33 


2.89 


1.89 


5.74 






60 


12.24 


1.17 


14.07 


2.82 


1.86 


6.06 






100 


15.04 


0.88 


12.09 


2.11 


1.60 


7.86 


E 


90 


1 


6.08 


0.38 


20.59 


0.92 


2.72 


1.85 






10 


7.82 


0.36 


19.69 


0.86 


2.60 


2.69 






20 


9.42 


0.33 


18.67 


0.80 


2.46 


3.51 






30 


10.75 


0.31 


17.84 


0.74 


2.36 


4.19 






50 


11.75 


0.28 


17.01 


0.68 


2.25 


4.85 






70 


13.07 


0.26 


16.00 


0.62 


2.11 


5.78 






110 


14.71 


0.23 


14.87 


0.55 


1.96 


6.97 


F 


100 


1 


8.74 


1.58 


22.41 


3.81 


2.96 


3.34 






30 


14.22 


1.34 


19.27 


3.34 


2.54 


6.64 






70 


17.62 


1.38 


16.67 


3.33 


2.20 


8.78 






100 


19.81 


1.28 


15.31 


3.08 


2.02 


10.17 






210 


24.43 


0.83 


12.50 


2.00 


1.65 


13.49 



Note. — Models A-F arc the same as in Table 2. Columns 4 to 9 are the following masses in 
units of Mq : (4) the remnant mass for the ease of O.O7M0 ^^*^Ni ejection {M^^-^/Mq), (5, 6) 
the total Mg and O masses in the ejecta, (7, 8) the '"^^^Ni masses when the ejecta has the solar 
abundance ratios for Mg/Fe and O/Fe (''^*^Ni[]yjgypj,]=o ^^i*^ ''^'^Ni[o/Fc]=o)> respectively, and (9) 
the upper limit to the ejected '''^Ni mass when the mass-cut is set just above the Fe-core shown 
in Table 2 (^^Niup) as a function of progenitor mass M and explosion energy £^51 = erg. 
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Log p_, (g cm 3) 

Fig. 1 . — Evolution of the central density and tem- 
perature for the 30 and 90Mq models. 
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Fig. 2. The enclosed mass (Af,.) - radius rela- 
tion of our pre-SN progenitors. The parameters 
and some properties of the models A-F are shown 
in Tabic 2. For each model the location of the top 
of the Fe-core is shown by a filled circle. i?Ni is 
the upper radius for the region where Ni^^ is dom- 
inantly produced as a function of the deposited 
energy E* given by the Equation (5) in the text. 




Fig. 3. — The enclosed mass (M^) - density (p) 
relation of our pre-SN progenitors. The masses of 
the models are 25, 30, 50, 80, and 100 M©, from 
bottom to top respectively. 
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Fig. 4. — The enclosed mass (M^) - electron mole 
fraction (Vg) relation of our pre-SN progenitors. 
The masses of the models are 25, 30, 50, 80, and 
100 Mq, from left to right respectively. Solid lines 
represnt 30 and SOM© models. 
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Fig. 5. — The post-explosion abundance distribu- 
tions for the 50Mq models (upper two panels) and 
IOOM0 models (lower two panels) when the ex- 
plosion energy is -E51 = 30 (left) and -E51 = 100 
(right.) 
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Fig. 7. — The upper hmit of ^^Ni masses produced 

by core-collapse SNc with various masses and ex- 
plosion energies. ^^Niup shown here are the ejected 
^^Ni mass when the mass-cut is located at the top 
of the Fc-corc. If mass-cut is larger the ejected 
mass will be smaller. 



Fig. 6. Same as Figures 5 but the horizontal axis 
is replaced by the radius of the progenitor, i?ini 
(km). The location of i?Ni calculated by the equa- 
tion (5) are also shown by large triangles. Here, we 
note that Rm appears to be a little over-estimated 
because the temporal energy absorption due to the 
matter dissociation is not taken into account in the 
equation (5). 
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